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The proposed effort involved development of numerical procedures for 
efficient solution of free vibration problems of spinning structures. 

In particular, the following R & D task has been completed as proposed. 

Task 1. An eigenproblem solution procedure, based on a Lanczos 
method employing complex arithmetic, has been success- 
fully developed involving 

(i) Formulation of numerical procedure 

(ii) Fortran coding of the algorithm 

(iii) Checking and debugging of software 

(iv) Implementation of the routine in the STARS program. 

Task 2. A graphics package for the E/S PS 300 as well as for the 
Tektronix terminals has been successfully generated and 
consists of the following special capabilities 

(a) A dynamic response plot for the stresses and 
displacements as functions of time 

(b) A menu driven command module enabling input of data • 

on an interactive basis. 
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Task 3. The STARS analysis capability has been further improved by 

implementing the dynamic response analysis package that provides 
information on nodal deformations and element stresses as a 
function of time. 

A number of test cases were run utilizing the currently developed algorithm 
implemented in the STARS program and such results indicate that the newly 
generated solution technique is significantly more efficient than other exist- 
ing similar procedures. 

The following section presents some details of the solution results. 
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NUMERICAL EXAMPLES 


The results of numerical analysis of two representative structures analyzed 
by the recently developed general purpose finite element computer program 

Q 

STARS , incorporating the current solution techniques are presented in this 
section. Thus, the first example relates to a square cantilever plate rotating 
with a uniform spin rate along an arbitrary axis. Figure 2 depicts the plate 
with a 6 x 6 finite element mesh, the edge along the X-axis being clamped and 
the structure having the following basic structural data: Young's modulus 

(E) = 10 7 . thickness (t) = 0.1, side length (l) = 10, Poisson's ration (y) = 
0.3, mass density (p ) = 0.259 x 10 . A free vibration analysis of the non- 

spinning structure was initially performed yielding a first natural frequency 
) value of 3.60 (expressed in non-dimensional parametric form defined in 
Table I) . Subsequent analyses of the structure were performed for the cases of 
= 0.8 mi as well as for a resulting spin vector Q R = 0.8 03 x having components 
Q x = ^y = = 0*8 u>i//3 T Numerical results for these analyses are presented 

in Table I. 

A coupled helicopter rotor-fuselage system^ freely floating in space, as 
shown in Figure 3, is chosen as the next example. Associated varying stiffness 
and mass distributions are suitably approximated for the discrete-element 
modeling of the structure. Free vibration analysis was performed with the 
rotor spinning at 10 rad/sec ( Qy = 10) . The results are presented in Table 
II along with the results for the corresponding non— spinning case. Associated 
mode shapes are shown in Figure A, that correspond to the non— spinning rotors. 


CONCLUDING REMARKS 


A generalized numerical formulation has been presented for the effective 
evaluation of nodal centrifugal forces in various finite elements due to any 
arbitrary spin rate. This, in turn, enables the derivation of element geo- 
metrical stiffness and centrifugal force matrices Kq and K°, respectively, 
which are vital for the free vibration analysis of spinning structures. The 
paper also presents the details of an improved eigenproblem solution procedure 
for the efficient free vibration analysis rotating structures, that effect 
considerable saving in solution effort when compared with the earlier procedure 
presented in Reference 1. A recently developed general purpose finite element 

Q 

computer program, STARS , incorporated current advances such as the techniques 
described herein. It has proven to be useful as an analysis tool for the 
solution of complex practical problems. 
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FIGURE CAPTIONS 


Figure 1 Triangular thin shell element spinning along an arbitrary axis 
Figure 2 A cantilever square plate spinning along a specified axis 
Figure 3 Coupled helicopter rotor-fuselage system 


Figure 4 


Helicopter mode shapes 


Table I. SPINNING CENTILEVER PLATE: natural frequencies for various spin rates 
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(b) Structural mass distribution 



Structural stiffness distribution 


Figure 3 
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